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ABSTRACT 

Several correlations among parameters derived from modelling the high-energy properties of GRBs 
have been reported. We show that well-known examples of these have common features indicative of 
strong contamination by selection effects. We focus here on the impact of detector threshold truncation 
on the spectral peak versus isotropic equivalent energy release (-Epk-^iso) relation, extended to a large 
sample of 218 Swift and 56 HETE-2 GRBs with and without measured redshift. The existence of 
faint Swift events missing from pre-Swift surveys calls into question inferences based on pre-Swift 
surveys which must be subject to complicated incompleteness effects. We demonstrate a generalized 
method for treating data truncation in correlation analyses and apply this method to Swift and pre- 
Swift data. Also, we show that the Ep^-E 7 ( "Ghirlanda" ) correlation is effectively independent of the 
GRB redshifts, which suggests its existence has little to do with intrinsic physics. We suggest that a 
physically-based correlation, manifest observationally, must show significantly reduced scatter in the 
rest frame relative to the observer frame and must not persist if the assumed redshifts are scattered. 
As with the E^-E ^ correlation, we find that the pre-Swift, bright GRB E p k-Ei so correlation of 
lAmati et al.l (|2006h does not rigorously satisfy these conditions. 

Subject headings: gamma rays: bursts — methods: statistical — cosmology: observations 



1. INTRODUCTION 

Correlations are pervasive in astronomy and gen- 
erally lead theory in allowing us to discover causal 
relationships between observed quantities. In the 
study of Gamma-ray Bursts (GRBs), multiple pow- 
erlaw relations among GRB observables have been 
used to unc over the intrinsic physics of GRBs them- 
selve s (e.g., lEichler fc L cvinson 2004; Yamaza ki et al.l 
2004 iRees fc Meszarosl 12005): iLevinson fc Eichierl 120051: 



Lamb, Donaghv. fc Graz iani 20051) and to use GRBs as 



probes to the dista nt Universe (e.g., iGhirlanda et al.l 
120051 : ISchacfer 2007|). Several authors have critically ex- 
amined the limitations of such relations: understand- 
ing how well they potentiall y constrain physics given 
their form and sca tter ( e.g., IFriedman fc Blooml [20051 : 
ISchaefer fc C ollazzi 2(Ml. 

uncovering possible ev olutio n 
with cosmic time (e.g.. lYonetoku et all [200l O [200l 
realiz i ng the commonality o f outliers dNakar fc Piranl 
[20051 : iBand fc Preecel [20051 : iKaneko et all 12006ft . and 
- most fundamentally — determining how the exis- 
tence and form of the relations vary once spurious cor- 
relation imparted by selection ef fects is treated (e.g., 
iLlovd. Petrosian. fe Mallozzill2000l ). 

Recently, iButler et al.l (120071 hereafter B07) r eport ev- 
idence from Swift satellite (jGehrels et al.ll2004j ) observa- 
tions of GRBs that several of the correlations exhibit a 
wide scatter and shift in normalization toward the Swift 
detection threshold, suggestive of an origin intimately 
connected to the detection limits of pre-Swift satellites. 
Because relations intrinsic to the physical processes un- 
derlying GRBs should not be instrument-dependent, a 
broader investigation into the data from pre-Swift satel- 
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lites is crucial for determining whether and/or how the 
relations can be trusted to potentially constrain the 
physics of GRBs or cosmology. 

Here, we extend our critique to include the full B07 
catalog of Swift GRBs with and without measured red- 
shifts z (for 218 total Swift GRBs) and also to include 
HETE-2 GRBs with and without measured z (56 GRBs). 
The uniform B07 catalog is novel for deriving bolometric 
fluences for all Swift GRBs (between GRBs 041220 and 
070509), without requiring tight error bars on the spec- 
troscopic fit parameters. The sample is therefore flux 
limited by the sensitivity of the Burst Alert Telescope 
(BAT; iBarthelmv et al.ll2005D: b ecause the X-ray Tele- 
scope (XRT; IBurrows et al.ll2b05| ) localizes to few arcsec- 
ond precision nearly all BAT GRB afterglows, additional 
flux limits associated with afterglow localization and host 
galaxy detection (important for example when consider- 
ing the GRBs with measured z) are likely not present or 
are far less important as compared to stronger flux limits 
imposed in pre-Swift surveys. 

A study simila r to th is has recently been conducted by 
IGhirlanda et alj ((2008) comparing a smaller sample of 
bright Swift GRBs to the Beppo-SAX s ample (although 
only t he sub-samples with measured z). [Ghirland a et all 
(2008) find that faint GRBs detected by Swift (and 
HETE-2; Section [2]) are missing from the Beppo-SAX 
sample, and there is no clear explanation for the miss- 
ing data in terms of either the SAX trigger threshold 
or a (likely higher) threshold resulting from a demand 
that tight error b ars be derived in the GRB spectral 
modelling. As in Ghirlan da et al.l (|2008h . we also fo- 
cus on the correlation between the isotropic equivalent 
energy j^j so and the peak in the yF v spectrum E ^ 
(|Llovd. Petrosian. fc Mallozzill2000l : lAmati et al.ll2002h . 

Selection effects in h igh-energy GRBs obser vables are 
well-documented fe.g.. [Lee fc Petrosian! [l996l. if rarely 
treated. For faint GRBs, an expected departure of the 
true distributions of flux, duration, etc., from the ob- 
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served distributions is expected given the steep GRB 
number density versus pea k photon flux F relation 
(dN/dF oc F- a , a ~ 2; e.g.. IPreece et al.ll2000D . which 
places most GRBs in a given sample near the detection 
limit — or near some other imposed sample cutoff at a 
higher flux level — with a narrow logarithmic dispersion 
of ln(10)/(l - a) ~ 0.4 dex. 

In Section [2] of this paper we study the correlation of 
Epk.ohs with the energy fluence S oc F E p ^ i, s dT roughly, 
where dT is the burst duration, for which it is clearly 
important to compensate for the expected non-intrinsic 
correlation of fluence with -Epk.obs that arises essentially 
because GRB detectors are photon counters and not 
bolometers. Similarly structured correlations are then 
discussed in Section [3] 

Two basic approaches have been attempted to com- 
pensate for GRB flux limits: (1) forward folding 
of a model GRB rate density, evolution, and lu- 
minosity function through a model of the detec- 
tor response (e.g.. iGraziani. Donaghv. fc Lamb! 120051 : 
lLamb. Donaghv. fc Grazianil l2005h . and (2) applying 
non-parametric statistical methods to the data and the 
observed flux limit for each data point. In Section [21 
we review and apply numerical methods which treat the 
observed data truncation directly according to (2) , while 
we plan to approach the B07 data via path (1) in a fu- 
ture paper. In order to potentially better understand 
the origin of pre-Swift correlations not strongly present 
in the Swift sample, we explore the significance of these 
correlations after rejecting GRBs from the Swift sample 
which would not have been detected by HETE-2 (Section 
0). 

Extending our critique to the iGhirlanda et alj (|2008l ) 
and similar pre-Swift samples of bright GRBs with red- 
shifts (e.g.. lAmati elaTl [200l [2009 : iSchaefel 12001 - 
where the data truncations are poorly known and the 
GRBs come from multiple instruments of varying sensi- 
tivity — we study (Section [3]) how the relations trans- 
form from the observer to the source frames. We show 
that redshift dependence in the correlations is generally 
weak, which is expected for a correlation that arises due 
to detection selection effects but also continues to have 
low scatter in the source frame. We present tests to un- 
cover whether the redshift dependence in other correla- 
tions is similarly weak, possibly providing circumstan- 
tial evidence that these correlations are non-intrinsic. 
This general approach, despite its limitations, can poten- 
tially find application in other population analyses where 
thresholds are important. 

2. FLUX LIMITS AND THE E PK -E ISO CORRELATION 

In Figure [TJ we reproduce the narrow (< 0.2 dex) 
scatter in E p k ~ E ®^ observed in the recent review by 
lAmati et al.l (120061) for data f rom multiple missions. Fol- 
lowing iNakar fc Piranl (|2005h in realizing that E^/E^ 
reaches a maximum for z w 3.83, we assume this redshift 
value for all GRBs without measured redshift to place the 
entire B07 sample on the plot (black points). Most (67%) 
best-fit Swif t value s are below the lower ler red-dotted 
lAmati et all (|2006l ) line. Nearly half (41%) are below 
the red dotted line at 90% confidence, indicating a clear 
preference for a lower flux normalization. We determine 
90% confidence error bars directly for E^/E p k, given 
the observed spectrum, rather than attempting to prop- 



agate errors on the covariant quantities Ei so and -E p k- 

If we plot the BAT threshold corresponding to 
the best-fit points, we can see that agreement (right 
side of the plot) is partly dictated by low-sensitivity, 
while many sensitive observations are strongly incon- 
sistent. This finding is in excel lent agreement with 
the BATSE studies o f faint GRBs (|Nakar fc Piranll2005l : 
iBand fc Pree^l2^^ .Jliere_ is no clea r separation in 
Figure Q] between t he lAmati et al.l (|2006D -consistent and 
lAmati et al.1 (|2006f )-inconsistent points, precluding obvi- 
ous (non-circular) cuts which could create the semblance 
of close consistency. 

In order to determine the most accurate threshold for 
each burst, we utilize the observed GRB spectrum and 
time profile directl y instead o f employing a model for 
the threshold (e.g., Band 2003) which would only utilize 
results from the spectral fitting (i.e., peak photo n flux, 
_E E k a obs, etc.) and possibly the burst duration (e.g., lBaiidl 
2006). After determining the energy band — taken to be 
15-350 keV for Swift but allowed to vary for HETE-2 (see 
Figure [5]) — and temporal region which maximizes the 
signal-to-noise ratio S/N, we divide the observed fluence 
by (S/N) /10. This assumes a background dominated 
light curve (i.e., one where J> 10 times more counts reach 
the detector from the diffuse sky background and poten- 
tially from non-imaged sources not of interest than from 
the source of interest), which is appropriate for all but the 
brightest handful of Swift GRBs, for which our threshold 
calculation takes into account the observed background 
counts. 

Here and below, we choose a cutoff S/N = 10, because 
this is where the observed S/N distributions appear to 
turn over. The turn-over indicates a drop in detection 
efficiency, because the number of fain t bursts is known 
from BATSE (e.g.. IPreece et al.|[2000h to increase at low 
flux levels. Because HETE-2 and Swift have intelligent 
trigger systems which seek to find the light curve region 
maximizing the S/N (e.g.. iBandl 120061) , the counts that 
define the S/N also define the fluence threshold, approx- 
imately. This threshold estimate is conservative in cases 
where the trigger system fails to find the optimal S/N 
region. 

Many of the BAT points in Figure [1] are likely lower 
than plotted, because most known GRBs have z < 3.83. 
Also, the lower limit error bars in the figure for most 
of the points are strongly influenced by a BATSE-based 
prior on -E p k,obs> which introduces a bias against the rare 
but not-unprecedented (e.g. from BATSE) inference that 
£ P k,obs > 300 keV (see, B07). 

Studies including some bright BAT GRBs (e.g., 
lAmati et al.1 120061 : iSchaeferl l2007f ) do not report the 
strong levels of inconsistency because the studies exclude 
the faint or hard GRBs for which tight -E p k,obs constraints 
cannot be obtained. However, for the simple question 
of consistency/inconsistency, tight -Ep^obs error bars are 
not required. 

It is interesting to note that short-duration GRBs tend 
to appear in the lower left of Figure [TJ populating a re- 
gion where GRBs of longer duration generally cannot 
be detected. Only one short-duration GRB is detected 
at a flux level consistent with the pre-Swift E^/Epk 
level. As for long-duration GRBs, short-duration GRBs 
tend to be detected near threshold. Appart from whether 
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In the early afterglow-era, these non-parametric 
methods were applied to cor relations between 
GRB sourc e frame quantities (iLlovd fc Petrosianl 
iLlovd. Petrosian, fc Mallozzil 120001: 
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Fig. 1.— The majority (> 67%) of Swift GRBs (black cir- 
cles) prefer a lo wer Ei BO -E D ^ normal ization relative to the pre- 
Swift GRBs in lAmati etaD ((2006:) (red circles). The Swift 
GRBs are sorted and indexed according to their 10cr limit- 
ing _E i ° S o/_E'pk values so that the GRBs detectable over the 
broadest range in the y coordinate are plotted toward the 
left of the plot. Swift GRBs with measured z are marked 
with a small yellow circle. For GRBs without measured z, 
we assume z = 3.83 , which maximizes E^/Ep^ (see, also, 
INakar fc Piranl 120051 ). Swift short GRBs (T 90 < 3 s) are 
marked with a cyan X. At these z, many (41%) of the Swift 
sample GRBs exhibit -E°n/-E P k values i nconsistent (at 90% 
confidence) with the (lAmati et all 2006) relation (red hori- 
zontal lines). The lAmati et al.l (|2006f ) consistent GRBs typi- 
cally lie to the right of the plot, where E^ /E p y_ sensitivity is 
weakest, indicating that these GRBs would be undetectable 
at significantly lower fluence levels. 

.Epk correlates with E- lso , because of the impulsive nature 
and resulting low fluence over peak flux ratio for short- 
duration GRBs, short-duration GRBs are more likely 
than long-duration GRBs to be outliers to an £? pk --Ei S o re- 
lation. We caution, therefore, against using low E^ / Ep^ 
values to classify (as opposed to, e.g., Amati 2008) inter- 
mediate duration events as short-duration (i.e., potential 
binary merger) events. 

2.1. From Swift BAT to HETE-2 

One straight-forward way to diagnose whether a flux 
limit is truly the origin of inconsistency between the BAT 
and pre-Swift samples is to compare the BAT data to 
data from satellites with very different detection sensi- 
tivities. If the corresponding surveys both extend to the 
respective limiting flux limits, then powerlaw fits for each 
sub-sample should have appropriately different normal- 
izations. (In pre-Swift datasets considered alone as in the 
next section, the relative detection sensitivities generally 
do not vary widely, and a different technique must be 
adopted.) In this sub-section, we compare the Swift and 
pre-Swift samples using a well-known approach which di- 
rectly treats the flux limits. 

The application and refinement of statistical methods 
in surv ival analys i s was pi oneered in BATS E GRB stud- 
ies by iPetrosianl (|1993l ). iLee fc Petrosianl (| 19961) . and 
others. By restricting comparison to "associated sets" 
of data - those GRBs detectable above the estimated 
thresholds of each other GRB in the set - it is possible 
to study distributions of and correlations among observ- 
ables independent of the thresholding, without making 
assumptions on the nature of the GRBs below threshold. 



iLlovd-Ronning fcTTl amircz-Ruiz 2002; Koccvski & Liand 
120061 ). The study of ILlovd. Petrosian. fc Malloziil (|2000l) . 
in particular, is important for (1) demonstrating that 
the correlation between isotropic equivalent 7-ray peak 
flux Ljso and -E p k likely arises as a result of the detection 
process while also (2) discovering a possibly intrinsic 
correlation between the isotropic equivale nt energy 
release E lso and E p ^ (also. lAmati et a l. 2002). 

The Kendall's tk statistic employed in these stud- 
ies reports the fraction of concordant (i.e., correlating) 
data minus the fraction of discordant data (i.e., anti- 
correlating) . Use of the statistic to rule on a correla- 
tion's strength derives from maximum likelihood prin- 
ciples, and the r^-test i s non-paramet ric and has max- 
imum statistical power (|Kendalll 119381 ). A simple, ele- 
gant, and rigorous ex tension of the Tj<--test acc ounts for 
data truncation (see, lEfron fc Petrosianl 119991 ): tabula- 
tion is restricted to the concordant or discordant data 
pairs which are detectable above (or below in the case 
truncations from above) each other's limits. In this way, 
both lower and upper limits on either the x or y variable 
can be treated, without the need to make assumptions 
regarding the missing data. Assumptions, sometimes dif- 
ficult to uncover and possibly also hard to justify, regard- 
ing the missin g data must generall y be made with other 
methods (e.g.. iGelman et al.ll2004T) . 

A major advantage of the truncated r^-test (which 
employs the limits directly) is that it is in general much 
easier to calculate a limit for an individual burst than it 
is to derive a general limit in E mo valid for all E p ^. In 
fact, the general limit can never be calculated precisely, 
because the flux limit (e.g., a satellite trigger threshold) 
is not a function of i?; so , but of the observed photon flux 
in some band pass over some shorter time integration. 

In the observer frame (Figure [2]A) , there is a highly 
significant apparent correlation between Sboi and -E p k b s 
for both HETE-2 (t k = 0.47, 5.2a) and Swift (t k = 
0.45, 9.3a); see Figure[3] Here, Sboi is the energy fluence 
calculated i n the source frame 1-10 4 keV band if z is 
known (e.g.. IAmati"eT al. 2002) or in the observer frame 
1-10 4 keV band if z is unknown. The significance of the 
correlation becomes modest (tk — 0.27, 4.4cr) when we 
account for the Swift threshold. 

That some significance remains indicates some amount 
of true correlation between -Epk.obs and .Sboi; although, 
this residual correlation is strongly affected by the flux 
limit. Quantifying the slope, scatter, etc., of this possi- 
ble residual correlation and its origin in local intrinsic or 
population-wide evolution effects will be the subject of 
a future study. Here, we are interested in the expected 
appearance of the correlation in the Swift survey rela- 
tive to pre-Swift surveys. The significance of the Swift 
(t k = 0.21, 2.4cr) or HETE-2 (j K = 0.16, 1.2cr) correla- 
tion becomes marginal if we account for the HETE-2 de- 
tection threshold (see, Figure^). To apply the HETE-2 
threshold to Swift data, we fit a second order polynomial 
to the HETE-2 threshold and evaluate that curve at the 
Swift E'pk.obs points. 

In the source frame and restricting to the 63 normal, 
long-duration GRBs with z in B07, we find tk = 0.59 
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(6.8<r) for the -E p k-£iso relation. However, considering 
only associated sets of observations above Swift thresh- 
old, the correlation strength drops strongly (tk = 0.28, 
2.4cr). The similarly strong decrease in correlation sig- 
nificance for the source and observer frame correlation 
can be understood by noting that the correlations involve 
very similar data: observer frame points above and below 
the HETE-2 threshold transform to source frame points 
with a very weak dependence on the GRB redshifts. The 
red dashed curve in Figure 03 — corresponding to the 
trajectory one event follows as its z is varied — illustrates 
how z only weakly affects where points fall on the plot 
(see also, Section[3|). Points in the extreme Bottom- Right 
of Figure or Top-Left of Figure 03 c annot be made 
consistent at any z (Figur e [TJ see also, iNakar fe Piranl 
l200l IBand fe Preecdl2005l) . 

Randomly sampling redshifts from those observed by 
Swift, we find that (80 ±3)% of Swift events above the 
HETE-2 threshold are consistent with the pre-Swift re- 
lation (solid and dotted lines). Only (20 ±3)% of Swift 
event below the HETE-2 threshold are consistent. Most 
(60%) of the events under the HETE-2 threshold are in- 
consistent at any z. We find similar results — (62 ±4)% 
consistency above HETE-2 threshold and (15 ±3)% con- 
sistency below HETE-2 threshold — if we assume a very 
different (toy) z distribution and draw ln(z) from a unit 
normal distribution. The best-fit E v \cE\ sa curve from 
B07 is a factor w 2 higher than the pre-Swift curve (solid 
and dotted lines in Figure 03). 

3. THE ANATOMY OF A CORRELATION 

We now turn our attention to a set of correlation tests 
useful for GRB samples with poorly understood flux lim- 
its. As mentioned above, pre-Swift afterglow-era satel- 
lites do not have typical sensit ivities that vary by more 
than a factor of a few (see also,|j3and 2003); we find that 
diagnosing the reliability of a correlation plotted on a log- 
log scale for these data cannot effectively be done by look- 
ing at variations in the normalizations alone. Moreover, 
the truncations tests discussed in the preceding section 
cannot easily be applied because the flux limits for each 
GRB in most pre-Swift surveys are unpublished. We seek 
to understand how the regression slope and scatter of the 
relation fits vary with assumptions on the nature of the 
relations. 

For simplicity, if we assume that all GRB fit param- 
eters have the same error, the errors can be ignored. 
Consider observer frame quantities X — log (Sboi) and 

Y = log (£pk,obs), for example. These translate to source 
frame quantities X' = X + bZ = log (£7i so ) and Y' = 

Y + aZ = log (-Epk), for example. Here, Z = log (1 + z) 
and a = 1, b m 3 (accurate to 20% for 0.5 < z < 5, 
representing the bulk of known GRB redshifts). 

If we assume an origin in the observer frame, which 
we hope later to rule out, cross terms involving 
the observables and z can be ignored, and the ex- 
pected source- frame linear regression slope (see, e.g. 
iBevington fc Robinson! [2003T ) is: 
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Here, (.) denotes an average after subtracting away the 
mean in X or 2, and m Q b s is the slope in the observer 
frame. Equation [T] gives the slope of the apparent source 
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Fig. 2. — From data truncation to spurious amplification of 
correlation strength in the observer and source frames. Typ- 
ical error bars are plotted for points in the Top-Right cor- 
ner. (A) There is a strong (> 9cr significant) correlation in 
the observer f rame between -Epk.obs and Sboi for Swift (B07) 
and HETE-2 (|Sakamoto et alj|2005r ). The jagged, gray line 
shows the observed detection S/N relative to S/N = 10 for 
Swift GRBs, and the solid black line is a fit to those val- 
ues. Similarly, using the larger S/N value from either FRE- 
GAT E or WXM from the HETE-2 catalog (|Vanderspek et al.1 
2008), we determine a 10er threshold curve for HETE-2 bursts. 
If we consider associated sets (see, e.g.. iPetrosianl fl993f) of 
data points above the HETE-2 threshold (e.g., dashed line 
regions), the correlation in the Swift data is low (2.4a). The 
significance is still only 2.7a if we lower the estimated HETE- 
2 threshold by a factor of 2, an appropriate amount to ac- 
count for error in our threshold estimate (e.g., jagged line). 
HETE-2 events with measured z's are marked with an X. (B) 
Randomly sampling z's from Swift bursts, we show one iter- 
ation placing Swift points from (A) onto the _E p k--Eiso plot. 
In 10 4 iterations, only (20 ± 3)% of Swift event below the 
HETE-2 thr e shold are consistent with the relation found by 
lAmati et all (|2002j ). 
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stituting into equations Q] and [2 it follows that X 
3/41og(S bol T jct ), Y = log(£ pk , obs ), a = 1 and b « 3/2. 
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Fig. 3. — The inferred -Epk.obs-Sboi correlation strength for 
Swift and HETE-2, as measured with Kendall's tk coeffi- 
cient, drops strongly as we correct for data lost to first the 
Swift threshold and then the (higher) HETE-2 Sbo\ threshold 
(Section [2]). 

frame relation given the observed data, whereas deter- 
mination of the true source frame relation requires de- 
tailed knowledge of the GRB rate density and luminos- 
ity function to impute the missing data. Although re- 
cent progress has been made on measuring these (e.g., 
iKistler et a.1.1 l2008h . we do not attempt here to recon- 
struct m bs [fisrc.true] • Moreover, for comparison with 
most GRB studies which simply assume m srCi true = rn STC 
(e.g.. lFirmani et aT1l2005l ). we only require Equation[T] 

Relatively independent of which survey we use, 
(Z 2 )/(X 2 ) w 0.1, and m src rs m b s - This is the typi- 
cal expected behavior for GRB correlations where X or 
Y are typically much more broadly distributed than Z. 
There are, however, e xample correlations (e.g., consider 
Willingal e et al.l i2007T) where the observables X and Y 
are arranged such that m ob s ~ 0, and the scatter in Z 
relative to that in X can be important for defining the 
chance m src according to Equation [1] 

We define the scatter a in Y about the best-fit regres- 
sion line to be the root-mean-square (RMS) deviation. 
For the best-fit regression line in Equation [I] the ex- 
pected scatter is 

a s 2 rc = ((Y - m SIC Xf) + (a- m SIC b) 2 (Z 2 ). (2) 

Covariance terms between X and Z or between Y and 
Z — which generally act to decrease the scatter in 
the source frame as compared to the observer frame — 
are dropped because we assume the correlation is non- 
intrinsic. To be explicit, rearrangement of an equation 
like E p k oc shows that the observables should vary 
with z, and this implies cross terms in Equation [2] The 
cross terms go to zero if the correlation is due to (z- 
independent) observational effects. 

The first term in Equation [2] is bound from the bottom 
by <7o bs ; equality is obtained for m src = m bs- The sec- 
ond term in Equation [2] implies an increase in a 2 ^ over 
°obs unless, coincident ally, a = m SIC b. Interestingly, this 
"coincidence" is satisfied nearly or exactly for many of 
the known GRB correlations. 

As an example, consider the IGhirland a et al.l {2004) 
correlation: E p y oc E®- 7 , where E 1 oc E iso 9 2 ct is the en- 
ergy release corrected for beaming into a jet of angle (9j Ct . 
For the case of a uniform density medium surrounding 



the GRB (|Sari et alJ ll99flh. 0? t oc E£' 4 (1 + z)- 3 / 4 7fj 4 , 
where T ]ct is the time in the observer frame where the 
jetting effects are expected to become apparent. Sub- 



Hence, a 



2 b, as also demon s trated graphically in 



Figure [4jA- The IGhirlanda et all (|2004l ) correlation is, 
therefore, effectively independent of the measured spec- 
troscopic redshifts. 

3.1. Example Numerical Tests 

For the lGhirlanda et"aT1 (|2004h correlation data in Fig- 
ure 21^, the RMS scatter is (7 obs = 0.17 dex for the ob- 
server frame variables and cr src = 0.16 dex for the vari- 
ables corrected for cosmological distance and redshift. 
This is a 0.9a insignificant decrease according to an F- 
test. Likewise, we find that the source frame data are cor- 
related strongly with Kendall's tk — 0.7. However, data 
generated with randomized redshifts (e.g., red points in 
Figure [U^) exhibit larger tk in most (76%) of the simu- 
lations. Somewhat more encouraging, we find a marginal 
2cr increase in significance in <T src ,sim = 0.18 ± 0.01 dex. 
However, the increase is apparently dominated by 1-2 
outlier events in a fraction of the simulations; if we in- 
stead employ the median absolute deviation about the 
median of the fit residuals as a robust measure of scat- 
ter, 52% of simulations exhibit lower scatter than for the 
observed data. 

For comparison with the previous sectio n and consid- 
ering t he 47 normal, long-duration bursts in lAmati et al.l 
(2006) for the E pk -E iso relation, we find a significant 
change in RMS scatter and slope only if we include two 
events at Ei SO < 2 x 10 51 erg. Using the outlier-resistant 
measures of correlation (tk or scatter estimated via the 
median absolute deviation about the median) we find 
that 10% of simulations yield a better correlation than 
the observed correlation, independent of whether the two 
XRFs are excluded. Hence, there is only very weak 
(1.5cr) evidence from our tests to favor an intrinsic ex- 
planation fo r the E u k-Ei SO correlation using the bright 
lAmati et alJ (|2006h data. Exclusion of a small numbe r 
of outliers is common practice (e.g.. lAmati et al.ll2006l) : 
however, it may be reasonable to criticize us for allowing 
the outlier events to be flagged during and not before the 
simulations. 

As an interesting side-note wh ich may further suggest 
a paucity of intrinsic physics in IGhirlanda et al.l (|2004l ) 
type relations, the balancing with z is characteristic of 
similar re lations found assu ming instead a wind-stratified 
medium (|Nava et al.ll2006f ) o r even allowing j9j e t t o vary 
arbitrarily with Tj et and E- lso (|Liang fc Zh ang 2005) . The 
latter relation has the smallest <7 src , which leads us to 
suspect that the only reason that a very large number 
of correlations do not abound in the literature is that 
we typically restrict to those than can be interpreted 
physically, or, short of this, we require the observables 
to evolve with z in a physically plausible fashion (e.g., 

Tjet,src = Tjet/[1 + z]). 

Without these constraints, it is possible to invent 
highly-significant but absurd relations (e.g., rcdshiftcd 
Swift trigger number versus redshifted burst duration) 
which make a mockery of all correlation studies. This ef- 
fect al one is likely not sufficie nt to generate the tightness 
of the IGhirlanda et al.l ( 2004h correlation , howe ver. 

The tightness of the IGhirlanda et al.l (|2004h correla- 
tion, if it is not intrinsic, must stem additionally from the 
effects outlined in the previous sections for the function- 
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Fig. 4. — iGhirlanda et al.l (|2004f ) find a tight correlation 
between the beaming-corrected GRB energy re l ease E -y and 
-Epk- We also plot recent points from ISchaeferl ([2007) . Sug- 
gesting a non-intrinsic origin similar to the correlations in 
Figure [2] it is not necessary to know the z of an event to 
acceptably place it on the best-fit curve (dotted line). The 
dashed red curve shows the trajectory with z, holding the 
other observables fixed. The tightness of the correlation is 
unaffected whether true or random redshifts are assumed for 
translating the observer frame to source frame quantities. 

ally s imilar E^-E^q c o rrelat ion. Excluding GRB 970508 
as in rGhirland a et aTl (|2004f h we find that the -E p k-£i SO 
relation has only a slight increase in scatter (0.18 dex) 
relative to the E^-E^ relation (0.16 dex), considering 
the same data (Figure [4]). In fact, because the inferred 
Tjet values are very narrowly distributed and E~ is pro- 
portional to E- mo to a power less than unity, a E p ] i -E 1 
relation is potentially always tighter than a E p ^-E lso re- 
lation, depending on how the scatter is measured. 

4. DISCUSSION 



The B07 Swift BAT and ISakamoto efaTI (|2005h 
HETE-2 GRB samples both exhibit a statistically sig- 
nificant correlation in the observer frame between S\, a \ 
and £pk,obs- What creates these correlations? The cor- 
relations follow the trigger threshold limits (Figure [2j\) 
in both cases, and the correlation significances drop pre- 
cipitously when we account for the flux limits (Section 
EJ. Therefore, flux limits must play a strong role shap- 
ing the observed correlation, perhaps giving them most 
of their statistical significance. 

What does this tell us about the source frame -Epk- 
Ei so relation? Because the redshift dependence in the 
transformation from £pk,obs-<Sboi to E p ^-Ei so is weak, the 
source frame correlation is likely to have the same origin 
as the observer frame relation (Section [3|). Indeed, we 
observe a strong decrease in the source frame correlation 
significance when we account for the Swift threshold. 

Finally, what can be learned about the E p ^-Ei so re- 
lation as it appears in Swift data relative to the pre- 
Swift E p k-Ei so relation? The normalization appears to 
be strongly instrument-dependent, and this suggests the 
pre-Swift normalization is defined largely by pre-Swift 
satellite flux limits. B07 find that the Swift E p k--Ei SO re- 
lation has a lower flux normalization and more scatter 



relative to pre-Swift Ep^-Eiso relations. We show above 
that we can raise this normalization (and decrease the 
scatter about the relation) to a level consistent with the 
pre-Swift E p \.-Ei so relation by imposing a heightened flux 
limit corresponding to detection by a satellite of HETE- 
2-like sensitivity. 

Strictly speaking, from these observations we can only 
rule out the existence of a narrow relation between E pk 
and i?i S o, while some physical correlation between the 
quantities may be present at high flux levels (see, e.g., 
the paucity of bright events in Figure [2(3). In any case, 
future satellites more sensitive than Swift are expected 
to further shift and broaden the E p ^-E lso relation into 
an inequality. 

More strongly, because we can probably attribute the 
slope (see, Section [3]), scatter, and normalization of the 
correlation to selection effects related to a photon flux 
F cutoff and the functional correlation b etween -E Dk . obs 
and S bo\ oc FE p \. t0 \, s dT (Section 1; also. iMassaro et aTl 
2007), an intrinsic explanation for the correlation in the 
Swift data may be unnecessary. Our results also call into 
question any intrinsic explanation for the correlation in 
pre-Swift data or in su rveys only including t he bright- 
est Swift events (e.g., Sakamo to et al.l [2008). because 
an E p k-Ei so relation should be instrument-independent. 
The existence of faint Swift events (in the observer frame) 
mostly missing from pre-Swift surveys shows that we 
have a poor understanding of the observational selec- 
tion effects which tru ncate the earlier data (see, also, 
IGhirlanda et~afl [20081 1. 

We could be incorrect in drawing these conclusions if: 
(a) our spectral fits are systematically incorrect and there 
are, in fact, no GRBs at flux levels below HETE-2 thresh- 
old levels, or (b) we have estimated the HETE-2 thresh- 
old incorrectly by a factor > 2. The accuracy of the 
spectral fits is addressed extensively in B07, where direct 
consistency is established relative to observations from 
Konus-WIN D or Suzaku of burs ts also detected by Swift. 
We note that lBellm et al.l ([2008D have independently ver- 
ified the statistical analysis in B07 for several bursts also 
observed by RHESSI. Assumptions in B07 regarding the 
GRB spectra (see, Section [2]), useful in compensating 
for the narrow BAT bandpass when determining bolo- 
metric flucnccs, are conservative in that they err toward 
large Ef^/Epk values. It also seems unlikely that the 
HETE-2 threshold should have any bearing on the cen- 
soring of Swift GRBs, considering also that many authors 
have reporte d BATSE GRBs detected below HETE-2 flu- 
ence limits dNakar & Pirad 12001 iBand fc Preecel 120051: 
iKaneko et al.ll2006h . 

For (b), we have shown that the Sboi-Epk.obs cor- 
relation significance, after accounting for the HETE-2 
threshold, does not increase if we lower our threshold es- 
timates by a factor of two. A larger error than this on 
our part for many GRBs is extremely unlikely given the 
straight-forward nature of the threshold calculation (Sec- 
tion H}. Moreover, such a substantial increase in HETE- 
2 sensitivity would make HETE-2 as sensitive as Swift, 
which is demonstrably incorrect given the dramatic dif- 
ference in the GRB localization rate of the two mis- 
sions (~20/yr for HETE-2 and ~90/yr for Swift). Given 
the mean peak photon flux to energy fluence ratio from 
ISakamoto et~aiT(|2005l 10 8 E pk ,ob s ph s" 1 erg" 1 keV), 
our HETE-2 threshold estimate is within 50% of that 
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estimated in lBandl (|2003D . 

Finally, there is strong indication that flux limits as- 
sociated with spectroscopic redshift determination are 
important to estimate and consider wh en restricting to 
GRBs with measured redshift (e.g, as inlGhirland a et aTl 
120081 ) . In Figure[2]we mark the HETE-2 events with mea- 
sured redshift; these are on average a factor two brighter 
relative to threshold than the events without measured 
redshifts. 

The quality of the HETE-2 GRB localization and pos- 
sibly also the brightness of the optical transient depend 
on the brightness of the GRB; both effects contribute to 
whether the afterglow and host galaxy can be detected. 
Additional complicated selection ef fects re l ated t o deter- 
mining the GRB z are discussed in lBlooml (|2003l) . These 
are expected to be less important for Swift due to arc- 
second X-ray localizations. 

We stress that a correlation is not de facto intrinsic 
simply because trigger thresholding can be ruled out as 
influencing the surveys. Other flux limits may domi- 
nate (and probably do dominate when the samples are 
restricted to GRBs with measured redshift). A crucial 
step in utilizing pre-Swift data to rule on the nature of 
correlations is to establish completeness in the surveys. 

4.1. New Correlations: Handle with Care 

To gauge the importance of systematic effects in this 
area of research, we have isolated from the literature two 
paths that likely have lead to apparently highly signifi- 
cant correlations: (1) selection effects truncate the data 
in various ways and the "missing" data are not treated; 
or (2) partial correlation with a hidden variable or vari- 
ables is ignored. In this paper, we have studied type (1) 
errors. However, both type (1) and (2) errors can be 
unmasked using the tools outlined above. 

Type (2) errors potentially arise in studi es which em- 
ploy one correlation (e.g., Lag-Luminositv. iNorris et all 
2000) to infer z for another correlation with similar vari- 
ables (e.g., i?pk-Luminosity) , without controlling for par- 
tial correlation with the variable in common (see, e.g., 
Llovd-Ronning fc Ramirez-R uiz 20021; iKocevski fc Liang) 
20061) . 

For example, from BATSE satellite observations where 
the detection threshold is very well characterized in terms 
of peak luminosity, we know that the -E p k-Luminosity 
correlation is largely formed in th e detection process 
(|Llovd. Petrosian. fc Mallozzil I2000D . However, this po - 
tential strong bias is ignored when I Yonetoku eTaD(2004) 
consider a purely intrinsic iJpk-Luminosity to derive a 
significant evolution in the rest-frame properties of GRBs 
(e.g., the correlation between Luminosity and z). The 
potential instrumental origin of the E^ k-Luminosity re- 
lation is also ignored bv lSchaeferl (|2007l ). 

We recommend testing potential new correlations (or 
known correlations not explicitly mentioned here) by 
proving: 

1. Increase in RMS scatter cr b s > c src is statistically 
significant according to an F-test, and 

2. the correlation scatter or significance determined 
using outlier resistant measures strongly decreases 
when the z's are randomized and the correlations 
are recalculated, and 



3. the observables, grouped to one side of the correla- 
tion equation, vary with z as predicted. 

These tests are to establish basic confidence in an in- 
trinsic nature for the correlations and are not merely to 
establis h that a correlation can be use d to estimate z 
re.g..lLl l200a ISchaefer fc Collazzi|[200l . 

For these tests to be accurate, however, it is also nec- 
essary to identify selection effects acting on the data. 
Flux or fluence values should be present down to the 
established survey completeness level. Well-established 
methods can then be applied to compensate for data 
truncation (Section ^ and to con trol for partial correla- 
tions (e.g., [Akritaii&ISiebir^lilll) . Covariance between 
the measured quantities, if present, must a lso be treated 
(e.g.. lLee fc Petrosian|[l99a ICabrera et al.ll2007l B07). 

5. CONCLUSIONS 

We show above for GRBs observed with multiple satel- 
lites that four example correlations reported in the lit- 
erature have features (instrument-dependent normaliza- 
tions, weak z-dependence, etc.) indicative of strong 
contamination by or even an origin in selection ef- 
fects. Contrarily, there is no widely-accepted, non-a- 
posterior explanation for the correlations in the source 
fram e. Also, the a-posteriori theoretical e xplanations 
(e.g jEichler fc Levinsonl2004lSchaeferl2007h fix the cor- 
relation normalizations by requiring GRBs to have one 
intrinsic spectrum or to be st andardizable candles (e.g., 
a nar row E 1 distribution; see, iBloom. Frail, fc Ku lkarni 
120031), a hypothesis no longer well-supported by the data 
(see. IKocevski fc Butlerll2008l) . 

The common z-independence of GRB correlations is 
either very odd or damning. As we discuss, this is one 
characteristic of a tight, apparent source- frame correla- 
tion which arises purely due to selection effects. For cor- 
relations between luminosity and some other measured 
quantity that does not depend on luminosity distance, 
z-independence is likely not a consequence of intrinsic 
physics, because the GRB cannot possibly know how the 
distance to the observer should vary with z. However, z- 
independence might be expected for correlations between 
luminosi ties (e.g., the recently r eviewed Ei so -Lx.w cor- 
relation, iNvsewander et al1l2008l ). 

We note that a requirement of z-independence trivially 
explains the relative slopes of the E p \,-E lso and Ep^-E^ 
relations (Section [3]), something theory can apparently 
do as well but only with reference to complicated bi- 
ases stemming from GRB beaming ([Levinson fc Eichlerl 
l200l. 

As an obvious point, we caution that the balancing 
of 1+0 terms on both sides of the correlation equations 
(from which the approximate z-independence arises) im- 
plies the correlations cannot be used to infer z. More 
speculatively — because redshift balancing allows a non- 
intrinsic observer-frame correlation to appear in the 
source frame as a low-scatter correlation — this balanc- 
ing may have played a role in the discovery of poten- 
tial intrinsic correlations. If the correlations currently 
known have been selected in place of correlations that do 
not balance redshift, and if the fitting of GRB properties 
and the pruning of outliers has been conducted with these 
correlations in mind, then using these correlations to con- 
struc t a Hubble diagram and test concorda nce cosmology 
(e.g.. iGhirlanda et al J 120051: lSchaeferll2007t ) is circular. 
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We stress that these redshift dependency problems — 
and our intrepetation of them — are unique to and po- 
tentially only characteristic of high-z objects like GRBs 
(as opposed to, e.g., SNe) where there is essentially no 
low-z calibration, and luminosity distance must generally 
be calculated using z and assuming a cosmology model. 

Small sample sizes may also have played an important 
role in allowing for tight apparent correlations through 
over-fitting. The correlation invol ving Luminosi t y, E p u , 
and T 45 duration for 22 GRBs in lFirmani et all i|2006). 
which does not appe ar to suffer from re dshift balanc- 
ing (see Figure 7 in iFirmani et al.l [2006( 1 . is an inter- 
esting case. When looking at larger datasets, however, 
there appears to be increased scatter and no statisti- 
cally significant improvement in scatter re lative to the 
. Epk-Luminosity ([Collazzi fc Schaeferl |2008[ ) or E v ^-E lso 
(|Rossi et"al"1l2008l) correlations. Investigation of the cor- 
relations in the largest possible datasets is, therefore, 
critical. 

The GRB community is no longer starved for data. 
The next critical step toward uncovering intrinsic cor- 
relations is to combine all available data by estab- 
lishing sample comple teness in pre-Swift surveys (e.g., 
iGhirlanda et al.1 [2008?) and treating the dominant flux 
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truncations using methods like those outlined above. 

In looking for new relations relevant to the physical 
processes underlying GRBs, and to avoid an inherent 
difficulty in deciding an origin for the correlations in the 
source or observer frames (Section [3]), it may be impor- 
tant to choose observables less-broadly distributed than 
the characteristic range in l+z or to abandon observables 
with strong and complicated truncations (i.e., fluxes 
or fluences) altogether. Functional correlations in l+z 
could be minimized by choosing observables (e.g., pow- 
erlaw indices) which do not vary explicitly with z. Most 
directly and circumventing all concerns raised above, we 
should establish calibration for GRBs at z <J 0.1. It is, 
therefore, also crucial to understand whether GRB prop- 
erties evolve with redshift. 
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